Are symmetric tidal streams possible with long-range dark-matter forces? 
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The unique dynamics of the tidal disruption of satellite galaxies is an extremely sensitive probe of 
long-range interactions between dark-matter particles. Dark-matter forces that are several percent 
the strength of gravity will lead to order unity changes in the ratio of the number of stars in the 
leading and trailing tidal streams of a satellite galaxy. The approximate symmetry of the stellar tidal 
streams of the Sagittarius dwarf galaxy would thus exclude attractive dark-matter forces greater 
than 10% the strength of gravity which would entirely eliminate the leading stream. However, recent 
simulations suggest that dark-matter forces 100% the strength of gravity could completely strip the 
stellar component of Sagittarius of its dark matter, allowing for the subsequent development of 
symmetric tidal streams. Here we argue that these simulations use inconsistent initial conditions 
corresponding to separate pure stellar and pure dark-matter satellites moving independently in the 
host galaxy's halo, rather than a single disrupting composite satellite as had been intended. A 
new simulation with different initial conditions, in particular a much more massive satellite galaxy, 
might demonstrate a scenario in which symmetric tidal streams develop in the presence of large 
dark-matter forces. This scenario must satisfy several highly restrictive criteria described in this 
paper. 
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I. INTRODUCTION 

Understanding the nature of dark matter (DM) re- 
mains one of the most important outstanding problems in 
cosmology. If DM is a new fundamental particle, it may 
experience new non-gravitational interactions that would 
manifest as a violation of the weak equivalence principle. 
Such DM forces should be observationally constrained on 
all scales where DM influences structure, from the sub- 
galactic out to the Hubble radius. In important early 
work, Frieman and Gradwohl p], Q recognized that DM 
forces would violate the equivalence principle, and con- 
sidered cosmological constraints on DM forces of the form 



(1) 



where G is Newton's constant, m is the mass of a DM par- 
ticle, r is the separation between particles, A is a "screen- 
ing length" beyond which the force is suppressed, and /3 
is a dimensionless parameter characterizing the relative 
strength of this force compared to gravity 35j. Obser- 
vations at the time of Frieman and Gradwohl's work ex- 
cluded DM forces with (3 > 1.3 and A > 0(1 Mpc). 

DM forces received renewed attention after it was re- 
alized that they were generically produced by a coupling 
to dark energy Q and could arise naturally in string the- 
ory 0, IE). They also might agree better with observa- 
tions than the standard cosmological constant plus cold 
dark matter ( ACDM) scenario of structure formation Q . 
Kesden and Kamionkowski (0,@|, hereafter KK06) ex- 
amined the effects of DM forces on tidally disrupting 
satellite galaxies, where the large difference between the 
mass m sa t of the satellite galaxy and that of the host 
galaxy Mh os t establishes a hierarchy in the relevant en- 
ergy scales of the problem. Assuming that the screening 
length A > d, an attractive DM force increases the satel- 



lite's orbital energy per unit mass 

Gil/host 
H/nrh — : 
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by a factor /3/DM,host/DM,sat; where d is the distance be- 
tween the satellite and host galaxies and /dm is the DM 
mass fraction in the satellite or host. A tidally disrupted 
star will end up in the leading (trailing) stream if its 
binding energy per unit mass is less (greater) than the 
typical energy 



E h 



(3) 



where r is the radius of the satellite galaxy. Equating 
the change in orbital energy induced by DM forces to 
the binding energy at the tidal radius 
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we find that DM forces will produce an order unity 
change in the ratio of the number of stars in the lead- 
ing to those in the trailing tidal streams when 
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KK06 applied their constraints to the Sagittarius (Sgr) 
dwarf galaxy, a Milky Way satellite approximately d ~ 17 
kpc from the Galactic center discovered in 1994 Q af- 
ter the publication of Frieman and Gradwohl's work. 
The Sgr dwarf has leading and trailing stellar tidal 
streams stretching hundreds of degrees across the sky 
that have been observed extensively by the Two-Micron 
All-Sky Survey (2MASS) [nj and Sloan Digital Sky 
Survey (SDSS) [ll(. The proposed mass-to-light ratio 
(M/L)s gr = 14 - 36(M/£) Q of the Sgr dwarf suggests 
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that it is DM dominated, as is the Milky Way at the 
perigalactic and apogalactic distances of the Sgr orbit, 
about 15 and 60 kpc respectively [HI]. The current bound 
mass of the Sgr dwarf is (2 — 5) x 1O 8 M0, implying 
Wsat/A^host — 10~ 3 . Inserting this into Eq. ([5]) along 
with /dm, host ~ /dm, sat ~ 1, we find that DM forces sev- 
eral percent the strength of gravity should induce order 
unity changes in the leading-to-trailing ratio of stars in 
the Sgr tidal streams. N-body simulations confirmed that 
(3 > 0.1 completely eliminated the leading tidal stream 
for a wide range of models of the satellite galaxy, host 
galaxy, and satellite orbit. KK06 thus conservatively ex- 
cluded (3 > 0.1, as the Sgr tidal streams are observed to 
be roughly symmetric. 

Keselman, Nusser, and Peebles ([HI, hereafter 
KNP09) have proposed an alternative scenario in which 
observations of the Sgr dwarf could be consistent with 
(3 ~ 1 DM forces. In this proposal, the differential ac- 
celeration experienced by stars and DM segregate them 
from each other well before the first simulated pericen- 
tric passage. The self-bound stellar remnant of Sgr is 
free of DM and thus immune to DM forces. It proceeds 
to develop symmetric tidal streams in the Galactic grav- 
itational potential well. The absence of a Sgr DM halo 
leaves the remaining stars more vulnerable to tidal heat- 
ing during pericentric passages. The resulting heightened 
velocity dispersion could be misinterpreted as evidence 
for DM if the virial theorem was misapplied to this highly 
nonvirialized system. KNP09 argues that just such a 
misinterpretation could explain the high observed mass- 
to-light ratio (M/L) SgT = 14 - 36(M/L) Q of Sgr, which 
would otherwise be incompatible with the complete seg- 
regation of stars and DM required by j3 ~ 1. 

Unfortunately, the initial conditions used in the simu- 
lations presented in KNP09 are not appropriate for de- 
scribing this alternative scenario for the formation of the 
Sgr tidal streams. These simulations begin with the satel- 
lite galaxy's stellar and DM distributions moving with 
the same velocity and centered about the same point in 
the host galaxy's halo. In the absence of DM forces, 
two distributions with the same location and velocity are 
guaranteed to be gravitationally bound to each other. 
However, DM forces change the orbital energy of the 
DM particles by an amount AE ~ /J/DM.host-E'orb that 
for (3 = 1 exceeds the energy E^ n binding these parti- 
cles to the satellite. KNP09 notes that DM forces induce 
"a nearly complete segregation of the stars and DM well 
before the first pericentric passage." For these initial 
conditions, the stars and DM are already widely sepa- 
rated in orbital energy before the simulations even begin, 
though they initially coincide in phase space. This point 
is dramatized by noting that if these simulations were 
run backwards in time, by symmetry the stars and DM 
would be segregated just as efficiently. These simulations 
actually describe the temporary coincidental convergence 
of two different satellites, one purely stellar and the other 
pure DM, moving independently on very different orbits 
in the host galaxy's halo. 



Proper initial conditions would show that this scenario 
cannot produce a purely stellar satellite with symmetric 
tidal streams like the Sgr dwarf. The segregation of stars 
and DM seen in these simulations will happen on the very 
first approach of the satellite galaxy to the Galactic cen- 
ter, long before dynamical friction can put the satellite 
onto an orbit resembling that of the current Sgr dwarf. 
The satellite will have fallen in from distances greater 
than the host galaxy's virial radius, and having been ac- 
celerated by a (3 = 1 DM force will be moving at a ve- 
locity well above the gravitational escape velocity of the 
system. Once segregated from their DM halo, the satel- 
lite stars will be moving on a hyperbolic orbit and will 
be ejected from the host galaxy after a single pericentric 
passage. 

We have identified three conditions that must be satis- 
fied for the KNP09 scenario for the formation of the Sgr 
tidal streams to be viable: 

1. The DM force must be strong enough that differ- 
ential acceleration in the Galactic halo fully segre- 
gates the stars and DM prior to the formation of 
the tidal streams. 

2. The stellar density in the satellite core must be high 
enough for it to remain a bound object after being 
pulled free from its DM halo. 

3. The DM force must not accelerate the satellite be- 
yond the gravitational escape velocity of the host 
galaxy before the satellite stars arc pulled from 
their DM halo. 

We will describe our models of the host and satellite 
galaxies in Sec. [TTJ We will then determine whether 
these models satisfy the three conditions given above in 
Sec. HIl We will find that while (3 = 1 DM forces are 
more than strong enough to fulfill condition [T] they are 
too strong to satisfy condition [3] Star formation is also 
unlikely to be efficient enough in dwarf galaxies to yield 
stellar densities that satisfy condition^ A summary and 
some concluding remarks will be provided in Sec. IIVI 

II. GALACTIC MODELS 

A. Host Galaxy 

As we shall see in Sec. IIII Al for /3 > 1 stars will be- 
come unbound from the satellite DM halo at galacto- 
centric distances d much greater than the scale lengths 
of the Galactic bulge or disk. For convenience, we will 
therefore model both these components as point parti- 
cles at the Galactic center of masses Mbuige and Mdisk 
respectively. Following KNP09, we will adopt values of 
M bu i ge = 3.4 x 10 10 M© and M bu i ge = 1O U M for these 
quantities. KNP09 uses the density profile 

PDM,host,KNp(rf) = U2 + fo2)2 ( 6 ) 
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for the Galactic DM halo, with i>haio = 131.5 km/s and 
b = 12 kpc. We will assume spherical symmetry for the 
both the host and satellite galaxies, implying that for 
each component X a mass 



M x (d) = 4tt f 
Jo 



r 2 px(r) dr 



(7) 



is contained within a sphere of radius d. The DM pro- 
file of Eq. ^ is unsuitable for our purposes because 
-^DM.host(d) does not asymptote to a finite value at large 
d, implying that an escape velocity 



(8) 



cannot be defined. We will instead use a Navarro-Frcnk- 
White (NFW) profile QJ, [TJ 



Pdm(gO 



Po 



(r/b)(l + r/by 



(9) 



with the same scale radius b = 12 kpc and the nor- 
malization po chosen to contain the same Galactic DM 



mass within the initial apocentric distance d a 



80 



kpc simulated in KNP09. This implies a virial mass 
M Vj host = 1-05 x 10 12 M Q , virial radius d v = 206 kpc, 
and concentration c v ,host = d v /b = 17.2. all very reason- 
able values for describing the Milky Way halo. 



B. Satellite Galaxy 

Following the conventions of KNP09, we will model the 
stellar component of the satellite galaxy with a modified 
Hubble profile JQ] 



P*,satM = Pi 1 + 



-3/2 



(10) 



where r c = 0.55 kpc and p\ is chosen so that so that 
a stellar mass M* sat (rt) = 3.0 x 10 s M Q is contained 
within the tidal radius rt = 1.67 kpc. Unlike KNP09 we 
will assume that the stellar density drops to zero outside 
r t as our criteria will be applied to the satellite core. 
The satellite galaxy's DM halo will also use the NFW 
profile of Eq. ©, but with a scale radius r s = 3.7 kpc, 
virial radius r v = 18.65 kpc, and virial mass M VjSat = 
1.5 x 10 9 M Q . We will assume for simplicity that the DM 
density vanishes outside r v , as the DM at large radii is 
largely irrelevant. This choice of parameters implies 



M v . sat 



tit, 



0.2. 



DM 



(11) 



and therefore that star formation is 100% efficient in 
satellite galaxies. This assumption is difficult to reconcile 
with theoretical arguments suggesting that feedback and 
reionization suppress star formation in low-mass halos 
[l7| . but we will consider it as a starting point. 



III. CONDITIONS FOR SYMMETRIC TIDAL 
STREAMS 

The satellite and host galaxies described in Sec.HTlmust 
satisfy three distinct constraints if their interaction is to 
produce a tidally disrupting system with symmetric tidal 
streams in the presence of a j3 = 1 DM force. We will 
identify and describe these three constraints in the three 
subsections below. 



A. DM segregation 

KK06 showed that DM-dominated satellite galaxies 
will develop asymmetric tidal streams in the presence 
of DM forces. KNP09 argues that this can be avoided if 
nearly all of the DM is removed from the satellite galaxy 
before its stellar tidal streams begin to form. As the satel- 
lite galaxy descends deeper into the potential well of the 
host galaxy, both stars and DM will be stripped away 
by conventional gravitational tidal forces. This occurs 
because parts of the satellite at different distances from 
the center of the host galaxy experience different grav- 
itational accelerations. The weak equivalence principle 
implies that at the same position, all objects experience 
the same gravitational acceleration. However, DM forces 
violate the equivalence principle and individual stars, in- 
dividual DM particles, and the partly stellar, partly DM 
satellite galaxy will experience three different accelera- 
tions towards the host galaxy's DM halo. In a reference 
frame freely falling with the satellite as a whole, the in- 
dividual stars will experience a relative acceleration 

G/3/DM,sat-^DM,host(rf) 



a*-host [d) 



d 2 



(12) 



while the DM particles experience a relative acceleration 

G/3(l — /DM,sat)AfDM,host(rf) 



ODM-host 



(d) 



d 2 



(13) 



In the absence of DM forces (/3 = 0), these relative accel- 
erations vanish as required by the equivalence principle. 
Furthermore, stars will experience no acceleration rela- 
tive to a purely stellar (/dm, sat = 0) satellite galaxy and 
DM particles will similarly experience no acceleration rel- 
ative a purely DM (/dm. sat = 1) satellite. 

Tidal forces, supplemented by these relative acceler- 
ations, will succeed in disrupting the satellite galaxy if 
they are stronger than the attractive forces that act to 
keep stars and DM bound to the satellite. For stars this 
restoring acceleration is simply the gravitational attrac- 
tion 



t(r) = 



,sat [' )\ 



(14) 



between the stars and the satellite, while for DM particles 
this gravitational restoring force is supplemented by the 
DM attraction to the satellite's DM halo 



aDM-sat(^) 



G[M*, sat (r) + (1 + /3)M DM ,sat(r)] 



(15) 
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FIG. 1: The accelerations a experienced by stars and DM par- 
ticles a distance r from the center of the satellite galaxy and d 
from the center of the host galaxy. Left-hand panel: The lower 
solid black curve shows the gravitational attraction a*_ sat be- 
tween stars and the satellite galaxy. The short-dashed red and 
long-dashed blue curves show the contributions of stars and 
DM respectively to this gravitational attraction. The upper 
solid black curve shows the acceleration aoM-sat between DM 
particles and the satellite; the DM contribution has been dou- 
bled relative to a*_ sat for 1 -+- ft = 2. Right-hand panel: The 
short-dashed blue curve shows the tidal acceleration at the 
edge of the satellite's DM halo, the long-dashed blue curve 
shows the relative acceleration aDM-host, and the solid blue 
curve shows the sum of these two accelerations. The three 
red curves show the corresponding quantities for the stars. 
The satellite is fully stripped of DM when it first reaches a 
distance d scg — 32.7 kpc from the center of the host galaxy 
as shown by the vertical dotted green line in the right-hand 
panel. At this distance, the total relative acceleration experi- 
enced by DM particles (solid blue curve in right-hand panel) 
equals the maximum of the restoring acceleration (upper solid 
black curve in left-hand panel) which is located near the satel- 
lite's stellar core radius r c = 0.55 kpc (vertical dotted green 
line in left-hand panel). Tidal disruption of the stars does 
not begin until dud — 15.7 kpc (vertical dotted purple line in 
right-hand panel) when the total relative acceleration expe- 
rienced by stars (solid red curve in right-hand panel) equals 
the restoring acceleration for stars at the stellar tidal radius 
r t — 1.67 kpc (vertical dotted purple line in left-hand panel). 



These accelerations are shown in Fig. [T] for the galactic 
models described in Sec. [TI] with a ft = 1 DM force. The 
many curves in this figure are described in its long cap- 
tion. Disruption of the satellite's extended DM halo has 
already begun even at the host galaxy's virial radius d v = 
206 kpc. This can be seen by observing that the upper 
solid black curve aDM-sat in the left panel intersects the 
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FIG. 2: Upper panel: The radius r of satellite's DM halo 
(long-dashed blue curve) and stellar distribution (short- 
dashed red curve) as the satellite gets within a distance d 
of the host galaxy. Lower panel: The mass M of the satel- 
lite galaxy (solid black curve), its DM halo (long-dashed blue 
curve), and stellar distribution (short-dashed red curve) at a 
distance d from the host galaxy. The vertical dotted green 
curve shows the distance d seg = 32.7 kpc at which the stars 
and DM are fully segregated, while the vertical dotted pur- 
ple curve shows the distance d t id = 15.7 kpc at which tidal 
disruption of the stellar component begins. 



right boundary below where the solid blue curve aoM-sat 
in the right panel intersects the right boundary. As the 
satellite falls inwards, the radius r and mass M of its DM 
halo decrease as shown by the long-dashed blue curves 
in Fig. [21 Although the tidal forces initially dominate 
the disruption, the relative acceleration aDM-host exceeds 
the tidal acceleration for d < 147.6 kpc where the long- 
dashed blue curve crosses the short-dashed blue curve in 
the right panel of Fig. [TJ Disruption of the satellite's 
DM halo continues until d scg = 32.7 kpc, when the total 
relative acceleration of the DM exceeds the maximum of 
the restoring acceleration aoM-sat- At this distance all 
the remaining DM mass MDM,sat(?*c) = 1.6 x 10 7 M Q in- 
terior to the stellar core radius r c = 0.55 kpc where this 
maximum occurs is disrupted en masse. This occurs be- 
cause aDM-sat (f) is not a monotonically decreasing func- 
tion of r as a result of the cored stellar profile /0* )Sa t(r) 
of Eq. pTJ|) . The sharp drop in MrjM.sat leads to corre- 
sponding sharp drops in /dm, sat and a*-host in Eq. (Tl2"j) . 
explaining the cuspy nature of the total relative acceler- 
ation for stars shown by the solid red curve in the right 
panel of Fig. [1] This cored stellar profile is very conducive 
to the scenario of KNP09, as a cuspy profile would have 
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had a greater central restoring acceleration ODM-sat(O) 
and might have retained DM until stellar tidal disrup- 
tion commenced. The satellite galaxy, now entirely free 
of its DM halo, remains undisturbed until gravitational 
tidal forces on their own can disrupt it in the standard 
manner. This occurs at a distance dtid = 15-7 kpc from 
the host galaxy where the tidal forces exceed the restor- 
ing acceleration a*_ sa t(r t ) at the satellite's tidal radius 
r t = 1.67 kpc. 

Simulation S2 of KNP09 was designed to illustrate 
how symmetric tidal streams resembling those of the Sgr 
dwarf could come about even in the presence of = 1 
DM forces. It begins with the satellite at an apocentric 
distance of d ap = 84.8 kpc. Although the pericentric 
distance is not provided explicitly, if the streams are to 
resemble those of Sgr it must be close to Sgr's estimated 
pericentric distance of c? po ~ 14 kpc [l2|]. For this choice 
of an initial orbit we have 

d pc < dtid < d scg < d ap (16) 

which demonstrates one of the concerns with this simu- 
lation. Although the satellite begins the simulation with 
its DM halo intact out to the virial radius r v , this ex- 
tended halo should have been disrupted at much greater 
distances from the host galaxy. Since d pc < d sog , the 
satellite should not have any DM halo at all. These ini- 
tial conditions represent a temporary association of stars 
and DM moving independently in the Galactic halo, like 
ships passing in the night. To properly simulate this 
scenario, the satellite's orbit must evolve from that on 
which the DM was fully segregated (d pc = G? scg ) to the 
current orbit needed to produce the observed stellar tidal 
streams. 

Dynamical friction drives this orbital evolution, but 
was neglected in KNP09 since the streams develop 
on a few orbital times t or b and according to (l6l | the 
dynamical-friction time 

iDF ~ (^-] t orb ~ 10 3 i orb (17) 

is much longer than this for the Sgr-Milky Way system. 
In response to this paper, KNP09 performed a new sim- 
ulation with a "live" DM halo for the host galaxy that 
allowed for dynamical friction and a more massive satel- 
lite ?7i sa t /Mhost — 0.1 that would reduce the dynamical- 
friction time ioF- The satellite also began on an orbit 
with a greater apocentric distance d ap = 250 kpc to allow 
for orbital evolution. A detailed analysis of this new sim- 
ulation is beyond the scope of this paper, but it should be 
examined whether this live halo can adequately describe 
the dynamical friction given that the final satellite mass 
is m sat < 4 x 10 8 M Q and the individual N-body particles 
of the host galaxy's DM halo have masses of 10 8 M Q . 

The segregation distance d scg = 32.7 kpc for the 
Sgr dwarf is comparable to the current distances of the 
newly discovered ultrafaint SDSS dwarf galaxies Segue 
1, Ursa Major II, Willman 1, and Coma, all of which 



lie within 50 kpc of the Galactic center. These new ul- 
trafaint dwarf galaxies have a wide range of luminosities 
(10 2 < L/Lq < 10 7 ), yet share a common central den- 
sity M(r = 0.3 kpc) ~ 10 7 M@ [lj]. Simulations show 
that tidal heating cannot inflate the velocity dispersions 
of these satellites enough to eliminate the need for mas- 
sive DM halos [l|| . The satellites less luminous than Sgr 
will thus have higher DM fractions /dm. sat and corre- 
spondingly higher relative accelerations a*_host for stars 
according to Eq. (Ti"2"|) . The common central density how- 
ever implies that the restoring acceleration a*~ S at for the 
stars will remain largely independent of satellite lumi- 
nosity. The stars should therefore be fully stripped from 
their DM halos at greater distances than d seg = 32.7 kpc 
for the more massive Sgr dwarf. This is not observed for 
Segue 1, the closest of the new dwarf satellites, whose 
galactocentric distance of 28 kpc [2(| is well within its 
predicted d seg for (3 = 1. Although the ultrafaint dwarf 
satellites appear to place even tighter constraints on DM 
forces, if we restrict our attention to Sgr we see that a 
(3=1 DM force is strong enough to segregate stars and 
DM in Sgr at distances d seg > rftid before the stellar tidal 
streams begin to form. The first of our three criteria 
for Sgr to form symmetric tidal streams as proposed in 
KNP09 is thus satisfied. 



B. Intact Stellar Core 

Even if DM forces are strong enough to fully segregate 
a satellite galaxy's stars from their DM halo, the segre- 
gated stars may not survive as a self-bound object after 
being pulled free. The survival of an intact stellar core is 
essential to the scenario proposed in KNP09, as it is this 
remnant that subsequently develops the symmetric tidal 
streams observed today. The stellar core is threatened 
with destruction at three distinct stages: 

1 . The satellite galaxy may be stripped of all its stars 
before the stars and DM are fully segregated. 

2. The stellar component of the satellite may never 
have been self-bound, but merely held together by 
the gravity of the satellite's DM halo. 

3. Even if the stars are initially self-bound, they may 
be tidally disrupted by the satellite's DM halo as 
it is segregated. 

We will explore each of these threats in turn in the re- 
mainder of this subsection. 



1. Competitive Segregation 

We saw in the previous subsection that stars will be 
stripped from the satellite when the sum of the gravi- 
tational tidal acceleration and the relative acceleration 
a*-host exceed the restoring acceleration a*_ sa t, while 
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DM particles will be stripped when the same inequality 
is reached between aDM-host an d aDM-sat- This estab- 
lishes a competition between stars and DM: if the DM 
profile is sufficiently cuspy or the satellite is sufficiently 
DM-dominated, there may not be a stellar remnant left 
when the satellite finally reaches d se g. Returning to the 
accelerations plotted in Fig. [H we see that at large d the 
total relative acceleration for stars shown by the solid red 
curve in the right-hand panel was above the total rela- 
tive acceleration for DM shown by the solid blue curve. If 
the satellite's DM halo had been more massive or less ex- 
tended compared to the stellar component, /dm, sat would 
have remained close to unity and the total relative accel- 
eration for stars could have exceeded the maximum of 
the lower black curve in the left-hand panel before this 
happened for the DM at d sog . In this case, all of the stars 
exterior to this maximum located at r ~ r c would have 
been disrupted into asymmetric tidal streams. A small 
stellar remnant would remain because all of the stars in- 
terior to r c would have been stripped en masse as was 
the case for DM in our model of the Sgr dwarf in the pre- 
vious subsection. If the restoring acceleration a*_ sa t had 
been monotonically decreasing no bound stellar remnant 
would survive. 

Though we have estimates of the total mass and mass- 
to-light ratio of Sgr, the stellar and DM density pro- 
files themselves are very uncertain. It is worth exploring 
more generally which profiles will allow an intact stel- 
lar core to be segregated from its DM halo. A general 
DM density profile will scale as pdm(?") oc r~ a at small 
r; the NFW profile of Eq. © has a = I. The pro- 
file a = 1 is special in that the resulting acceleration 
Qdm(^) = GA/dmMA 2 oc r x ~ a will be independent of r. 
DM halos with a > 1 yield accelerations that diverge at 
r = and monotonically decrease with r. Density pro- 
files with a < 1 will produce accelerations that increase 
with r, before ultimately turning over at larger r since 
the halo must have finite mass. 

Of these three possibilities, the a = 1 profile is the 
most conducive to leaving a bound stellar remnant be- 
hind. DM halos with a > 1 will hold onto their innermost 
stars so tightly that d scg = will always be less than d t id 
at which stellar tidal disruption begins. Simulation S3 of 
KNP09, though still employing an NFW profile for the 
satellite halo, has an enhanced central concentration that 
mimics this possibility. When this simulation is run, the 
majority of stars are indeed stripped from the satellite 
halo while the inner stellar core remains bound. Satel- 
lites with a < 1 will have accelerations that vanish at 
r = 0, so even small differential accelerations adif will 
displace the stars from the center of their DM halo. As 
the satellite falls inwards, a*_host increases and the stars 
march up the DM potential well to ever larger accelera- 
tions a*_ sa t before breaking free from the satellite when 
they reach the global maximum. However, because the 
satellite stellar distribution has a finite spatial extent, its 
outer edges will spill over the global maximum before the 
center arrives there. This may cause additional stars to 



be stripped from the satellite before the eventual bound 
stellar remnant itself breaks free. 

If a was truly unity for an NFW profile, not just for 
r <C b, the acceleration &dm if) depicted by the long- 
dashed blue curve in the left-hand panel of Fig. [T] would 
be flat. If star formation was much less efficient as seems 
to be the case for less massive dwarf satellites, the DM 
contribution would dominate in a*_ sat and the lower solid 
black curve would be flat as well. In this case, the total 
relative acceleration would rise above a*_ sa t(r) simulta- 
neously over a wide range in r. The stars would slide 
out of their DM halo like "an egg sliding off a frying 
pan" [2l|, and the stars would be more likely to remain 
self-bound. In the competition to remain bound to the 
satellite galaxy, stars are more likely to defeat DM if stel- 
lar densities are high and centrally concentrated, while 
the DM density profile's logarithmic slope scales as a = 1 
over a wide range. 



2. Self-boundedness 

After being segregated from its DM halo, the stellar 
remnant may no longer have enough self-gravity to re- 
main bound. Unbound stars in the host galaxy's halo 
would develop into an orphan stream lacking a clear pro- 
genitor. One such system, apply dubbed the "Orphan 
Stream", was discovered serendipitously by SDSS [22| . 
This inspired simulation S4 in KNP09 which demon- 
strates just how readily such systems can be produced 
following the segregation of stars and DM by large DM 
forces. The possibility of making an orphan stream is a 
problem for efforts to explain the Sgr dwarf, which has 
a well defined bound core identified by a cusp in the M- 
giant density distribution spatially coincident with the 
globular cluster M54 This bound core serves as 

a marker to distinguish the leading and trailing tidal 
streams of Sgr, essential to the proposal of KK06 to con- 
strain DM forces by comparing the stars in each group. 

Whether the satellite stars remain bound after their 
DM halo is removed depends on their distribution func- 
tion fx.s&tir, w). If the stars are segregated from the DM 
on a dynamical time that is shorter than the relaxation 
time on which their individual energies change, one can 
use an impulse approximation to determine which stars 
remain bound. Without the negative contribution of the 
DM halo's gravitational potential, some portions of phase 
space that were previously bound will now have positive 
energies. The integral of /* lS at(^*, v) over this portion of 
phase space will determine the initial mass loss. With- 
out these stars the satellite's potential well will become 
even shallower, and the stars will rearrange themselves 
with further mass loss likely. Unfortunately, many differ- 
ent distribution functions will yield the density profile of 
Eq. p0[) . Distinguishing between them requires knowl- 
edge of the anisotropy of the satellite's velocity disper- 
sion, which is even more poorly constrained observation- 
ally than the density profile itself. 
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FIG. 3: Gravitational potentials $ as a function of the dis- 
tance r from the center of the satellite galaxy. The short- 
dashed red curve is the potential $» from stars, the long- 
dashed blue curve is the potential $DM<r t from DM interior 
to the tidal radius rt, and the solid black curve is their sum. 
The dot-dashed black curve is their difference $,-$dm whose 
stellar-mass- weighted average, appearing in Eq. 1|20[) . deter- 
mines whether the satellite stars remain bound after the DM 
halo has been removed. 

One can crudely estimate whether the satellite stars 
remain bound by applying the virial theorem to the initial 
system to determine its kinetic energy 

K = ~Wi (18a) 

= ~\j P*( r )[*.(»0+*DM<r t (r)] d 3 r. (18b) 

Here Wi is the initial gravitational potential energy of the 
stars, <3>*(r) is the gravitational potential of the satel- 
lite stars, and 3>DM<r t (0 is the gravitational potential 
sourced by DM located interior to r t . Spherical symme- 
try implies that only the mass interior to r t can accelerate 
the stars and thus contribute to their kinetic energy K. 
If the DM halo is removed quickly enough, this kinetic 
energy K will remain unchanged while the final potential 
energy will be 

W f = ^[p*{r)$*(r)d 3 ? (19) 

now that the DM halo has been removed. The total final 
energy Ef can be found by summing Eqs. (fT5)l and (fT^|) 

E f = K + Wf (20a) 
= \j P.(r)[*.(r)-*DM<r t (r)] d 3 f. (20b) 



If Ef > 0, the satellite as a whole will be unbound, 
though a bound remnant may still remain if the unbound 
stars carry away enough energy. Figure [3] shows $* and 
$DM<r, for the galactic models described in Sec. |TTJ The 
dot-dashed black curve shows that their difference, ap- 
pearing in the integrand of Eq. (|20|) . is negative for r <r t . 
Since p*(r) is positive definite, Ef must be negative im- 
plying that the stellar remnant after DM segregation re- 
mains self-bound. If star formation had been less effi- 
cient or the DM profile was cuspier, we could have had 
> $DM<r f for some r which might have driven Ef 
positive. Only star formation efficiencies approaching the 
theoretical maximum of Eq. (|11[) allow the possibility of 
leaving a bound stellar remnant behind after the stars 
are segregated from their DM halo. 



3. Tidal Disruption 

Even if the stellar remnant of the satellite is self- 
bound when initially displaced from its DM halo, it may 
be tidally disrupted by the satellite DM halo's gravity 
before it escapes beyond the satellite virial radius r v . 
This process is in some ways the reverse of traditional 
tidal disruption, in that the stars begin at the center 
r = of the satellite DM halo and are pulled to larger 
r by the relative acceleration, rather than beginning at 
a large distance d from the host and falling inwards to- 
wards d = 0. The tidal interactions between comparable- 
mass systems like the satellite stellar and DM distribu- 
tions are extremely difficult to model analytically, ne- 
cessitating the N-body simulations used by both KK06 
and KNP09. The best we can hope to do without such 
simulations is to compare the restoring acceleration at 
the edge of the stellar remnant, given by the stellar con- 
tribution to a»_ sat (r t ), to the tidal acceleration sourced 
by the satellite's DM halo given by the DM contribution 
to |a*_ sat (r) — a*_ sat (r + r t )\. The left-hand panel of 
Fig.[T]shows that the stellar contribution to a»_ sat (short- 
dashed red curve) in the range < r < r t is greater than 
the DM contribution (long-dashed blue curve) over the 
entire range < r < oo. This implies that the satellite 
DM halo cannot disrupt the stellar remnant during the 
segregation of the two components. As for the other two 
constraints considered in this section, this might not have 
been true for a cuspier DM profile or less massive stellar 
component of the satellite galaxy. 



C. Satellite Bound to Host Galaxy 

The final condition that must be satisfied for the 
KNP09 scenario for the formation of symmetric tidal 
streams is that the satellite stars remain bound to the 
host galaxy after they are segregated from their DM halo. 
This is a significant concern, as even in the absence of DM 
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FIG. 4: Velocities v of infalling stars, DM particles, and the 
satellite galaxy as a function of their distance d from the host 
galaxy. These velocities are for radial infalls beginning at 
rest from turnaround distances d ta equal to the host galaxy's 
virial radius d v = 206 kpc and 1.5d v . The short-dashed red 
curves correspond to freely falling stars, the long-dashed blue 
curves to DM particles, and the dot-dashed purple curves to 
the satellite galaxy which is partly stellar and partly DM. The 
solid black curve shows the gravitational escape velocity v eS c. 
The dotted vertical orange line shows the host virial radius 
d v , while the dotted vertical green line shows the segregation 
distance d sce . 
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gained during tidal disruption can sometimes unbind 
stars in the trailing tidal stream from the host galaxy. 
This energy E"tid is responsible for raising extended tidal 
tails in interacting galaxies, and has even been proposed 
as an explanation for the large fraction of hypervelocity 
stars with common travel times clustered in the direc- 
tion of the constellation Leo [23|. The satellite galaxies 
in simulations SI and S2 of KNP09 begin ~ 80 kpc from 
the host galaxy with velocities of 80 km/s, much less 
than the escape velocity v csc = 333 km/s at this dis- 
tance. Dynamical friction over several orbits is assumed 
to have placed the satellite on this orbit, similar to that 
of the Sgr dwarf, before segregation occurs and the tidal 
streams begin to develop. However, we saw in Sec. IIII Al 
that the stars will be segregated from their DM halo the 
first time the satellite gets with a distance d scg of the 
host galaxy. The satellite will be traveling at a greater 
velocity on this first approach and therefore its stars are 



more likely to become unbound from the host galaxy. 

The secondary infall of subsequent material onto an 
already collapsed structure has been studied extensively 
in a cosmological context 0, HH, HH. At early times 
this material recedes with the Hubble expansion, before 
reaching a maximum turnaround distance rfta an d falling 
back onto its host galaxy. Models of secondary infall can 
be used to predict the orbits of satellite galaxies [27j . 
and comparison with observations suggest d ta — 1 Mpc 
for the Milky Way [28:]. Surviving satellites like the Sgr 
dwarf are biased towards large turnaround distances dm 
compared to disrupted satellites that form the Galactic 
stellar halo [29[. 

In Fig. 0] we show the velocities acquired by stars, DM 
particles, and the satellite as they freely fall towards the 
host galaxy on radial orbits from the turnaround dis- 
tance. These estimated velocities are conservative, as 
the tangential velocities must be added in quadrature to 
determine the true total velocity. The velocities are cal- 
culated by assuming conservation of energy 



1 



aff(x) dx 



(22) 



where the different free-fall accelerations ag for the stars, 
DM particles, and satellite are 
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G[Af*,hoBt + (1 + /3)M DM ,ho S t(rf)] 
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(23a) 
(23b) 



a*-g(d) 
ar>M-s{d) 
a sa t-g(d) 



The DM fraction /dm. sat in Eq. (|23cp is calculated self- 
consistcntly as the ratio of the satellite's DM and total 
masses, given respectively by the long-dashed blue and 
solid black curves in the bottom panel of Fig. [2l Note 
that the difference |a sa t-ff ~a*-ff| of these free-fall accel- 
erations gives the relative acceleration a*_h os t appearing 
in Eq. (|12p . The difference |aDM-ff ~ Osat-ff| similarly 
gives the relative acceleration aDM-host in Eq. (|13[) . 

The attractive DM force proportional to (3 in 
Eqs. (|23bp and (|23c|) accelerates the DM particles and 
satellite to higher velocities than the stars can attain by 
gravity alone. This can be seen in Fig. 2] where the DM 
velocities Wdm-s given by the long-dashed blue curves 
and the satellite velocities u sa t-ff given by the dot-dashed 
purple curves exceed the stellar velocities f*_ff given by 
the short-dashed red curves. The satellite is DM dom- 
inated (/dm, sat — 1) at large distances, implying that 
i)DM-ff and fsat-ff are initially close together. As the 
satellite falls inwards and loses its DM halo, it becomes 
dominated by stars (/dm, sat — 0)and u 8a t-ff approaches 
i>*_fr . Crucially however, w sa t-ff remains above for 
the models of Sec.|H]it is even above the gravitational es- 
cape velocity v esc at d seg . The satellite is accelerated by 
the DM force on its DM halo while falling towards d seg , 
but once this halo has been removed the DM force cannot 
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help keep the stellar remnant bound to the host galaxy 
after the first pericenter passage. It should therefore be 
ejected from the host galaxy, not evolve onto the tightly 
bound orbit needed to produce stellar tidal streams re- 
sembling those of the Sgr dwarf as predicted by KNP09. 

Dynamical friction, described by the Chandrasekhar 
formula dHU 

' m satPhost - / n ,\ 
— 2 Vsat ' <^ 4 J 

dt v Lt 

might help keep the stellar remnant bound by reducing 
its velocity below the escape velocity. However Eq. (TTT1) , 
readily derived from the Chandrasekhar formula, sug- 
gests that dynamical friction will not be significant on 
an orbital time for a binary system with a mass ratio 
as small as that between the Sgr dwarf and Milky Way. 
Increasing the mass of the satellite m sat can increase the 
effect of dynamical friction, but introduces problems as 
well. A higher satellite mass will increase the tidal energy 
£"tid in Eq. (|2ip , which on its own can unbind tidal debris 
even in the absence of DM forces. The restoring acceler- 
ation aDM-sat will also increase with the satellite mass, 
reducing d scg and increasing the distance over which the 
DM force has the chance to accelerate the satellite be- 
yond the escape velocity according to Eq. (|2"2")l . An NT- 
body simulation will be required to determine which of 
these effects dominates for a particular model. 

Just such a simulation was performed in the latest ver- 
sion of KNP09 in response to the concerns expressed 
in this paper. This simulation has precisely the fea- 
tures necessary to avoid ejecting the stellar remnant 
of Sgr from the Galactic halo. The extraordinarily 
massive Sgr dwarf (M* iBat = 4 x 10 10 M Q , M DM , S at = 
2 x 10 11 M Q ) experiences strong dynamical friction ac- 
cording to Eq. (|24|) . and its large tangential velocity 
I'sat = Vcsc{d = 250 kpc) = 185 km/s allows time for 
this dynamical friction to remove orbital energy before 
the first close approach to the Galactic center. This sim- 
ulation demonstrates a possible scenario in which the Sgr 
dwarf could arrive on its present orbit in the presence of 
f3 ~ 1 DM forces, however several questions remain. The 
Sgr dwarf loses more than 99% of its stellar mass during 
the course of the simulation. Is this large a contribu- 
tion to the Galactic stellar halo consistent with obser- 
vations? Although the Galactic DM halo has now been 
made of "live" particles to allow for dynamical friction, 
the Galactic disk is still modeled by a static potential. 
Can the Galactic disk (M disk = (4.5 ± 0.5) x 1O 1O M 
[IB]) survive an encounter with such a massive satellite, 
particularly in the presence of DM forces? The simulated 
Galactic DM halo consists of 2 x 10 4 particles, each of 
mass 10 8 M Q . Is this sufficient to resolve the dynamical 
friction on the Sgr stellar core, which also has a mass 
A^*,sat — 10 8 M Q after losing 99% of its mass to tidal dis- 
ruption? While this new simulation is intriguing, further 
simulations and comparisons to observations are needed 
to confirm the validity of the KNP09 scenario. 

An alternative possibility to avoid ejecting the stel- 



lar remnant of the satellite is if it was accreted by 
the host galaxy as part of a larger group of satellites. 
Many of the Milky Way's brightest satellites appear to 
lie in a disk, as would naturally occur if these satellites 
were initially members of a group that was subsequently 
tidally disrupted in the Galactic halo [3l|. Lake and 
D'Onghia [32| specifically proposed that the Magellanic 
clouds formed the core of such a group that also included 
seven other Milky Way satellite galaxies including the 
Sgr dwarf. However, others have argued that this hy- 
pothesized group would have to be more tightly bound 
than observed dwarf-galaxy associations, and that inde- 
pendent accretion from a common filamentary structure 
could also explain the observed disk of satellites [33|. If 
the Sgr dwarf was accreted as the more massive mem- 
ber of a satellite-galaxy binary, it could lose energy to its 
lighter companion as the binary was tidally disrupted 
in the Galactic halo (34|. Such a "cosmic menage a 
trois" might conceivably leave the Sgr dwarf with a small 
enough orbital energy that its stellar remnant might re- 
main bound to the Milky Way despite the additional ac- 
celeration due to = 1 DM forces. 

IV. DISCUSSION 

KNP09 performed a series of simulations that sug- 
gested that the symmetric tidal streams of the Sgr dwarf 
galaxy might be consistent with DM forces comparable 
in strength to gravity, in contradiction with the claims of 
KK06. These simulations used static potentials for the 
Galactic bulge, disk, and DM halo. Such static potentials 
save computational resources and avoid artificial heating 
of the tidal streams, making them essential to attempts 
like those in [ljj to compare the observed and simulated 
velocity dispersions in the tidal streams. However, static 
potentials do not allow for dynamical friction and there- 
fore the satellite galaxy must begin on its current orbit to 
reproduce the observed tidal streams. Such a cheat was 
acceptable in [T^] and KK06, but fails in the scenario of 
KNP09 where the segregation of stars and DM should 
occur on the satellite's first approach within a distance 
dscg of the host galaxy. This first approach will likely 
occur on a very different orbit from the current one with 
apocentric distance e? ap ~ 80 kpc needed to produce the 
observed Sgr tidal streams. The initial conditions used in 
the first series of simulations in KNP09, with the satellite 
stars and DM coincident at d ap with the same velocity, 
actually correspond to already unbound stellar and DM 
distributions moving on widely separated orbits in the 
Galactic halo that happened to converge at t = when 
the simulations begin. 

We have considered the first approach of a satellite 
galaxy to the Galactic center in the presence of large 
DM forces, and have identified three preconditions for 
the creation of a purely stellar satellite that can go on to 
form symmetric tidal streams: 

1. The DM force must be strong enough that differ- 
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cntial acceleration in the Galactic halo fully segre- 
gates the stars and DM prior to the formation of 
the stellar tidal streams. 

2. The stellar density in the satellite core must be high 
enough for it to remain a bound object after being 
pulled free from its DM halo. 

3. The DM force must not accelerate the satellite be- 
yond the gravitational escape velocity of the host 
galaxy before the satellite stars are pulled from 
their DM halo. 

The extremely large mass-to-light ratios observed even in 
the cores of Milky Way satellites are difficult to reconcile 
with condition 2, and such efficient star formation con- 
tradicts our theoretical understanding of the vulnerabil- 
ity of these satellites to rcionization and stellar feedback. 
DM forces strong enough to satisfy condition 1 are likely 
to be too strong to satisfy condition 3 unless a further 
epicycle is added to give the Sgr dwarf an initial orbital 
energy much less than that expected for a newly accreted 
satellite galaxy. While the scenario proposed in KNP09 
for symmetric Sgr tidal streams is perhaps possible for 
certain initial orbits and galactic models, it requires that 
several restrictive assumptions be satisfied. Further sim- 



ulations are needed that self-consistently capture both 
the segregation of stars and DM and the subsequent de- 
velopment of tidal streams. 

Future observations of other Milky Way satellites will 
also severely constrain the existence of (3 ~ 1 DM forces. 
Once their proper motions have been measured, the peri- 
centers of their orbits can be determined. If these peri- 
centers are too close to the Galactic center, the DM forces 
they will experience during pericentric passages will be 
inconsistent with retaining their DM halos in seeming 
contradiction to the high observed mass-to- light ratios. 
More sophisticated simulations and observations by up- 
coming astrometry missions like SIM Lite thus provide 
a clear path towards closing remaining loopholes in con- 
straints on DM forces. 
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